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The dielectric a relaxation is the dielectric manifestation of the glass transition. In spite of 
this fact, the more commonly used models, the Arrhenius model and the Williams-Landel- 
Ferry model, do not take into account the structural state of the system to modelize this 
relaxation. In thermally stimulated discharge current (TSDC) experiments, the sample is out 
of equilibrium during most of the discharge ramp. Not surprisingly, the capability of these 
models to represent the data points at temperatures well below the glass transition temper- 
ature is very poor. To overcome this limitation, we have used the Tool-Narayanaswamy- 
Moynihan (TNM) model, that takes into account the structural state of the system, to 
modelize TSDC data. Although it is mostly applied to calorimctric and volumetric experi- 
ments we show how it can be employed on dielectric data. The numerical results support the 
applicability of the model and suggest how the dielectric and the structural relaxation times 
may be related. The TNM model turns out to be physically sound since the modelization 
of dielectric data gives also a reasonable structural kinetics of the system. The qualitative 
differences between the TNM model and the equilibrium models for different TSDC experi- 
ments is discussed. Experimental data coming from TSDC experiments where the thermal 
history of the sample is changed do not provide additional evidence but nevertheless are 
compatible with the parameter values that were obtained in the fits to the TNM model. 
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I. INTRODUCTION 

A glass is a material that behaves mechanically like a solid but has a disordered structure. In 
fact, glasses can be considered as overcooled liquids that have acquired amorphous rigidity. The 
most common way of making a glass is by cooling the material from the melt fast enough to avoid 
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crystallization (although many polymers do not crystallize at any cooling rate). If a liquid is cooled 
in such a way, molecules will rearrange progressively more slowly. Eventually the structure of the 
material will not be able to reach the equilibrium conformation fast enough to follow the cooling 
rate. This falling out of equilibrium is called the glass transition Q]. From this point of view, 
the glass transition would not be a true phase transition because all physical magnitudes would 
change in a continuous way. It has been suggested that it exists an underlying phase transition 
Q] but anyway the glass transition occurs across a narrow range of temperatures and therefore we 
can define a glass transition temperature (T g ). Due to the kinetic nature of the glass transition, 
T g depends on the cooling rate, but this dependence turns out to be rather weak [3]. The study of 
the glassy state properties has received a great amount of attention Q and it is nowadays a ver 
active subject of research, both about systems close to T g |3( or clearly below this temperature 

In a wide sense, an electret is a sample of dielectric material that produces a nearly permanent 
external electric field. This field results from permanent ordering of molecular dipoles or displace- 
ment and trapping of free-charge. The polarization and the relaxation processes of an electret are 
intimately related to the properties of the material. Therefore, an study of these processes can 
provide useful information about the material. A technique in wide use to perform such studies is 
thermally stimulated depolarization currents (TSDC) 

A classical TSDC study can be described in the following way: the sample is heated up to a 
polarization temperature (T p ) and polarized between two electrodes for a polarization time (t p ) by 
a polarizing electric field (E p ). Then it is cooled down to a temperature, called storage temperature 
(T s ), low enough so that depolarization takes place at a very low rate. If conventional polarization 
is employed, the polarizing electric field is applied not only during t„but also during a large portion 
of the cooling ramp. Alternatively, in window polarization (WP) 6} the electric field is switched 
off at a point of the cooling ramp placed a few degrees below T p or even at the beginning of the 
cooling ramp. 

Once the sample is at temperature T s , it is kept at this temperature for a storage time (t s ) and 
then it is heated at a constant rate while it is shortcircuited through an electrometer. The obtained 
displacement current is recorded as a function of temperature. In the thermogram that results, the 
relaxation processes can be seen as current peaks that form the TSDC spectrum. Electret behavior 
and TSDC have been widely described in the literature HQ. 

Despite their experimental simplicity, TSDC spectra of conventionally polarized electrets are 
not easy to interpret. On one hand, different mechanisms can contribute to the thermogram, such 
as space charge relaxation and dielectric relaxation. On the other hand, these mechanisms may 



not be elementary (i.e. they can not be described in terms of a single relaxation time). 

It is possible to resolve a complex TSDC spectrum into elementary spectra using the WP 
technique ^,0,0]. For this reason, WP is a useful technique to study relaxation mechanisms. As 
a consequence of employing WP only part of the mechanism remains activated. If the thermal 
window is narrow enough, the part of the mechanism that is activated will behave approximately 
in an elementary way. These spectra are well described by a single relaxation time. 

One of the relaxations that can be detected by TSDC is the a relaxation, which is the dielectric 
manifestation of the glass transition. In TSDC data, the a peak has its maximum at the dynamic 
value of T g This value is the one at which really takes place the glass transition during the 
leating ramp and usually lies some degrees above the static value of T g . Unlike in other relaxations 
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the Arrhenius law 



r(T)=T exp(j|) (1) 

does not work properly when the system is close to the glass transition because cooperative phe- 
nomena are involved [l^]. The temperature dependence is better described by a Vogel-Tammann- 



Fulcher (VTF) type expression 



r(T) = r exp 



E w 



(2) 



[R(T - r^) 

where is the temperature at which molecules would "freeze" and the relaxation time of the sys- 
tem would become infinite. It can be shown that this model is equivalent to the Williams-Landel- 
Ferry equation. In either case, r(T) has no longer its initial physical meaning, a characteristic 
time that represents the return of the excited mechanism to equilibrium. This is also due to the 
cooperative character of the mechanism. For this reason the pre-exponential factor has values that 
are not easy to interpretate in mechanical terms. In spite of this fact, it is still a valid empirical 
parameter to characterize the relaxation [si lisl]. 

Both expressions can be used to modelize TSDC data. Nevertheless, an important point is often 
missed: in eqs. Handle t(T) represents the relaxation time in equilibrium at a temperature T. To 
record a TSDC scan of the a relaxation, the sample has to be cooled well below T g and then heated 
over this temperature in a short amount of time. As a result, the configuration of the molecules 
that form the sample is far from equilibrium during most of the experiment. The departure from 
equilibrium depends on the kinetics of the system and the design of the experiment. 

There are a number of ways to extend these equations to describe the non-equilibrium glassy 



state 
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Many of them are phenomenological theories that have been often tested with 



calorimetric and dilatometric data [if]. \v\. Other models follow a more first-principles approach 
Occasionally, models that represent glassy behavior have been used to modelize TSDC 
data. For example, the coupling model j^J has been used for this purpose 21]. 

Polymethyl methacrylate (PMMA) is a polar, glassy and transparent polymer with multiple 
industrial applications. Its most common application is as a replacement for window glass because 
of its optimal chemical, mechanical and optical properties. Since this material has been intensively 
studied and the relaxations it can undergo have a rather simple behavior, it is also appropriate 
as a test field for new techniques such as space charge distribution measurements P] or second 
harmonic generation in guest-host systems |5j. 

The aim of this work is to demonstrate how the Tool-Narayanaswamy-Moynihan (TNM) model 
can be used to modelize data from the a relaxation of PMMA obtained by TSDC. Since this model 
takes into account the structure of the system, this approach highlights the relationship between the 
dielectric a relaxation, the structural relaxation and the glass transition. The physical implications 
and the coherence of the model will be discussed in order to evaluate its applicability to dielectric 
phenomena. 



II. EXPERIMENTAL 



The experimental setup for TSDC measurements is composed of a custom made convection oven 
placed inside a cold bath at —20 °C. The oven is controlled by an Eurotherm 2416 temperature 
programmer. Inside, there is a measurement cell with two electrodes that are applied to the 
sample in plane-parallel configuration. The current was measured by a Keithley 6512 electrometer. 
Temperature and current data were collected by an A/D converter card and processed afterwards. 

PMMA samples were cut from a commercial sheet of 1.5 mm thickness. The thickness of the 
samples was reduced to 0.75 mm by mechanical procedures. Circular aluminum electrodes of 1 cm 
of diameter were deposited in vacuum on both sides of the sample. Smaller samples were also cut 
for differential scanning calorimetry (DSC) experiments. 

DSC experiments were performed to determine the glass transition temperature at several heat- 
ing rates. We used a Mettler TC 11 thermoanalyser equipped with a Mettler-20 differential scan- 
ning calorimeter module. DSC curves were obtained from 15 mg samples sealed in aluminum pans. 
The heating rate of the scans was 20 °C/min, 10 °C/min, 5 °C/min and 2.5 °C/min. From a 
these results it can be extrapolated that T g ~ 110 °C/min. The dynamic value at 2.5 °C/min is 
T g « 115 °C/min. 



The TSDC experiments carried out follow the scheme represented in fig. ^ The values taken 
by the parameters shown in fig. ^ ar e listed in tabled All the ramps have a rate of 2.5 °C/min. 
Experiments start well above T g in order to clear the influence from any previous thermal history 
and to begin the experiment with the sample at equilibrium state. The sample is annealed for an 
annealing time (t Q ) at an annealing temperature (T a ). 

Next, the sample is polarized. In all the presented experiments T p = T a . Most of them use 
T p = 95 °C. This polarization temperature was chosen in order to isolate as much as possible the 
a peak from the space charge (p) peak, but without being too far from the optimal polarization 
temperature (T po ) of the a relaxation. In previous experiments with the same sample the T po was 
found to be 107.5 °C. A polarization potential of 800 V and a polarization time of 900 s (15 min) 
were employed in all the experiments. 

Once the polarization time is over, the field is switched off and a cooling ramp begins. This 
ramp continues until T s is attained. After t s the sample is short-circuited through an electrometer 
and heated again while the thermally stimulated discharge current is recorded. 



III. RESULTS AND DISCUSSION 



The TSDC thermogram that we will use to illustrate the application of the TNM model is shown 
in fig. [2 and has been obtained with the experiment #1 described in tabled The peak corresponds 
to the dielectric a relaxation. To analyze these data we will assume the Debye model. According 
to this model, the relaxation process presents an ideal-viscous rotational friction characterized by 
a single relaxation time (r) 



dP 1 , 



This magnitude is stro ngly dependent on the temperature, so we will refer to it as t(T). 

The Bucci method |22| is the framework in which the evaluation of the experimental value of 
t(T) has been obtained. This method is based in the fact that the Debye model given by eq. [3 is 
based in first-order kinetics. From this fact and other general considerations the expression 

T (T) = — — / I{T')dT' (4) 



vI(T) 

can be deduced. In this equation v is the constant heating rate and I(T) is the intensity of the 
recorded TSDC current. Given r(T) it is possible to recover the intensity (up to an arbitrary 



factor) from the relaxation time, using the formula 

/(T)ocex pj-I^-^-InKD]}. (5) 

The characteristics of this magnitude are best appreciated in a ln(r) vs. 1/T plot. To visualize 
better the correspondence between the features in the intensity plot and in the relaxation time 
plot, the intensity has also been plotted in fig. El using a reciprocal temperature scale. 

These data has been fitted to several models by \ 2 minimization. The routines used are de- 
scribed by other authors [23(. A point that deserves some attention is for which magnitude the 
maximum likeness should be achieved. A perfect model should give the same results independently 
of the magnitude used for the minimization, but this is obviously not our case. Although there 
are many possibilities, we have decided to perform our fits for two of them: / (as in |2lJ, |2J]) and 
ln(r) (as in [3]). The parameters resulting from the I fit, aside from giving better results for this 
magnitude, also represent better the behavior of the points placed near the maximum. Instead, 
parameters coming from ln(r) fits represent better the overall behavior of the experimental data. 
Anyway, the difference between the plots that have been computed with both sets of parameters 
gives a qualitative measurement of how appropriate is the model. This idea will be used throughout 
to assess the models. 

In many works (see for example |26|) the data fitted corresponds to the temperature range 
where I(T) is greater than I(T max )/2. In fig. [2] that would mean fitting the data points that lay 
between T{ = 95.5 °C and T2 = 110.5 °C. In this region ln[r(l/T)] has almost a linear behavior. 
Moreover, it has a slight concave curvature, so it is also possible to apply the WLF model. 

In fig. OK we compare the fits to an Arrhenius model with I data. The intensity plot that has 
been fitted to I data (point-dash-dashed line) represents fairly well the data placed between T[ 
and T2. On the other hand, the plot fitted to hi(r) data (point-dashed line) is very close to the 
previous one. The same happens in fig. Ob that represents ln(-r) data and plots obtained with the 
same parameters as in fig. EH- The plots obtained fitting to m(r) data and to / data are very 
similar. The parameter values obtained from these fits can be seen in table HTl 

In fig. HI we illustrate how the WLF model is even more successful than the Arrhenius model in 
this range of temperatures. Both kinds of fits, to I data and to ln(r) data, give almost the same 
parameter values (see table ITTTjl . The representation of both solutions is almost the same, either in 
fig. EH or in fig. EJ). 

We can see, then, that the equilibrium models, Arrhenius and WLF, represent finely the be- 
havior of the system in a narrow temperature range around the dynamic T g . Definitively, they 



are appropriate for this purpose, with the WLF even better than the Arrhenius model in this 
case. Let's now go on to consider what happens when a wider range of temperatures is considered. 
Instead of taking T[ as the lower limit, we will fit the data from T\ = 70 °C. In this way we include 
the data points that form the foot of the a peak but we still avoid the contribution that could arise 
from the (3 relaxation. It is not possible to enlarge the temperature range on the high temperature 
end, that we will keep at T2 = 110.5 °C. The reasons are that the accuracy of the integral in eq. 0] 
becomes too low and also that the contribution from space charge relaxation turns out to be too 
noticeable in the results. 

We will point out, in first place, that the WLF model can no longer be applied. There is an 
inflexion point near T[ in such a way that most of the ln[r(l/T)] curve has now convex curvature 
and the WLF model would give unphysical results, with > T g . Fig. El shows the results for the 
Arrhenius model (see numerical values in table [H]). It is clear that fits to / data and to ln(r) data 
give plots that differ much more than in the previous case. This is a clear sign that in this range 
of temperatures the model is no longer appropriate. 

Results obtained fitting ln(r) data do not represent well neither the relaxation time nor the 
intensity, as seen in fig. |31 The plots obtained fitting / data placed between T\ and T2 are very 
similar to the ones obtained fitting data points between T[ and T2 because, as mentioned earlier, 
this procedure takes more into account the points that are placed closer to the peak maximum. 
Nevertheless, in fig. OK it can be seen that the values of the intensity that lay between T\ and T[ 
cannot be modelized properly. 

A classical explanation would be that, even when WP is employed, several parts of the mecha- 
nism are polarized, giving rise to a distribution of relaxation times. According to this explanation, 
the thermogram can not really be treated as if it were elementary. 

In fact, it is true that the relaxation is distributed to some extent. Some experimental results 
can not be explained otherwise. Nevertheless, there are reasons to think that, even in this case, 
the Arrhenius model does not provide an appropriate basis to express the TSDC intensity. A 
broad distribution of relaxation times would be needed to modelize the region between T\ and T[. 
Instead, the data points between T[ and T2 can be modelized quite well using a single relaxation 
time. This suggests a narrow distribution, which is more likely in WP experiments. There would 
be no contradiction if the intensity at the beginning of the a peak corresponding to an elementary 
relaxation would rise up not as sudd enly as in the Arrhenius model. Otherwise, usual distribution 
functions such as Fuoss-Kirkwood j^Tj or Havriliak-Negami will not be able to give good 
results. Therefore, the problem remains more or less the same as before. 



There is a chance that this failure is due to the out-of-equilibrium state of the material during 
the TSDC ramp. In the following lines we will discuss how the TNM model can be applied to 
provide a more realistic profile of the TSDC intensity. To lighten up the discussion we will keep 
the formalism at the elementary relaxation level although it should be emphasized that the TNM 
model is not incompatible with the existence of a distribution of relaxation times. 

The TNM [3| phenomenological model has been used mostly to modelize calorimetric and 
dilatometric data. Within this model, the relaxation time depends not only on the present temper- 
ature but also on the thermal history of the sample. We will denote this relaxation time as r s and 
we will not suppose yet that it is identical or otherwise related to the dielectric relaxation time r. 

In the TNM model the structure dependence is taken into account through the fictive temper- 
ature (TV). The fictive temperature of a non-equilibrium system is equal to the temperature of an 
equilibrium system that has the same structure. The separation between temperature and fictive 
temperature dependence is introduced through a non-linearity parameter x (0 < x < 1). The 
structural relaxation time is given by 

ln(r.) = Info) + + (6) 

where E* is the apparent activation energy. The TNM model is reduced to the Arrhenius model 
given by eq. ^ either when the system is at equilibrium or in the limit when the non-linearity 
parameter is equal to 1. 

The kinetics of TV is determined by the first-order equation 



dT f _ T f -T 



(7) 



dt r s (T,Tf) 

In other words, we are assuming the same type of kinetics for both the dielectric and the structural 
relaxation. 

Eqs. 13 and are coupled. To solve the system, we define a reduced time 

xE* (l-x)E* 



z(t) = / exp 
Jo 



RT(t>) RT f {t') 
In terms of the reduced time, the solution of eqs. Eland [7| reads 



dt'. (8) 



T f {t)-T(t) = 0[z(t)][7)(O)-T(O)] 

-£^P-4>[z(t)-z(t)]dtr. (9) 

We have written the solution also in terms of the response function of the system 

</>(t) = exp (-—) . (10) 



This function determines the temporal behavior of the system in response to an external pertur- 
bation that has been applied such as an electrical field, a stress or a deformation. 

If T(0) > T g the first term of the RHS of eq. El vanishes because the Active temperature and the 
temperature are equal above T g . If the system is below T g at the beginning of the experiment, there 
is no easy way to know the initial fictive temperature. T(t) is a known function that reproduces 
the temperature from the beginning of the experiment to the end. It is possible to calculate t s (T) 
at any time of the experiment although only the t s (T) that corresponds to the TSDC ramp can 
be compared to the experimental one. 

Eqs. |H1 and El are as coupled as eqs. E] and and the integrals therein have to be evaluated 
numerically. To accomplish this, the time variable t is discretized and the integrals are calculated 
in the following iterative way. First, the fictive temperature of the previous time step is used to 
evaluate eq. El The resulting reduced time is introduced in eq.EJto calculate a closer approximation 
to the fictive temperature at the present time step. If the difference with the previous approximation 
to the fictive temperature is not small enough, the process is iterated. 

The time step taken has been of 5 s and Simpson's rule has been chosen to evaluate the integrals. 
Under these conditions, the convergence of the procedure that has been described is very fast and 
is achieved within few iterations, even if the precision required is of four significant digits. 

Up to this point, we know how to calculate the relaxation time of the structural relaxation (r s ) 
at any time for a given set of parameters (E*, t s q, and x) and a thermal history T(t). Now we 
should assume a relationship between r s and the relaxation time of the dielectric relaxation (r) . 

The simplest one is that both relaxation times are the same 

r(T)=r s (T), (11) 

but in our case it does not work. This can be seen in fig. El (numerical results in table llV|) . The 
model fits a little bit better to the experimental points because it has one more parameter and 
contains the Arrhenius model as a particular case but the same qualitative problems that we found 
for the Arrhenius model remain. In figEH> the parameters obtained fitting the TNM model intensity 
data give an almost identical curve to the Arrhenius model. Only the relaxation times close to 
the maximum are well modelized with these parameters (figEK)> just as happened using the same 
procedure with the Arrhenius model. The results obtained fitting the model to relaxation time 
data are as bad as in the Arrhenius model. The parameters obtained in this way are not able to 
reproduce the relaxation time data (figEJa), and obviously they perform even worse when they are 
used to reproduce the intensity (fig Eh)- To sum up, this hypothesis has more or less the same 



drawbacks as the Arrhenius model in spite of being much more involved. 

In terms of simplicity, the next relationship to consider would be a proportionality between the 
relaxation times 



t(T)<xt s (T). 



(12) 



This hypothesis implies the introduction of an additional parameter. This parameter can be 
expressed as another pre-exponential factor 



As within the previous hypothesis, the Arrhenius model is attained when x = 1, without the need 
of additional conditions. 

The results are shown in fig. |H1 (see table HVl for numerical values). The ability to match the 
experimental data is much better, either when we fit intensity data (continuous line of fig. EK) or 
relaxation time data (dashed line of fig. HjJs). In the case of the intensity, the model is able to 
modelize the points placed between T\ and T{, unlike previous models that could not cope with 
those points. The model can reproduce finely the position of the inflexion point and the curvature 
of each zone in relaxation time data. Another clear sign that the model performs better than the 
previous ones is that the curves obtained fitting / and ln(r) are very close, either in fig. HJ, or in 



Although this model is more accurate than the previous ones, some limitations can be noticed. 
The match of the theoretical curve with experimental data becomes worse as we consider points 
that are placed further below T g . This phenomena is reasonable since the TNM model is obtained 
by extending an equilibrium model (the Arrhenius model) and therefore it should be less valid as 
the system attains structural states further from equilibrium. The difference between the parameter 
values obtained fitting to intensity data and to relaxation time data is due to this disagreement. 

From the previous discussion it should be clear that the TNM model, together with the propor- 
tionality hypothesis, is able to fit TSDC experimental data with more accuracy than the Arrhenius 
and the WLF models. By itself this is not enough to consider that the TNM model provides an 
appropriate representation of the dielectric a relaxation. Better fits could be expected anyway due 
to the presence of more fitting parameters. On the following lines we will discuss the kinetics of 
the structure of the system within the TNM model to gain insight on its physical meaning and it 
will be checked that the kinetic behavior predicted by the model is the one that should be expected 
according to our current knowledge on physical aging. 



r(T) 



(13) 



TdO = TsO 



Ts(TY 



fig. ©3. 



In fig. we represent the evolution of the difference between Tf and T, for two different initial 
deviations. This magnitude can be considered as the deviation from structural equilibrium of the 
system. The time origin has been taken at the beginning of the cooling ramp that takes place after 
the first isotherm. 

The effect of the cooling ramp is to increase dramatically the deviation from equilibrium. The 
dotted lines have been introduced as a guide for the eye and represent slopes of 2.5 °C/min in the 
scale of the Tf — T plot. This is the rate of all the ramps in fig.0 From the plot follows that during 
the cooling ramp Tf tends to attain a constant value (about 90 °C for experiment #2, according 
to our calculations) and therefore the deviation will tend to grow at a rate equal to the one of the 
cooling ramp. 

The deviation does not decrease significantly during the storage time. For this reason during 
the TSDC ramp the deviation has its maximum value at the beginning of the ramp. Since the 
TNM is just a non-equilibrium extension of the Arrhenius model, this is the part of the ramp that 
will be harder to modelize. The dotted line allows us to see that, until t\, Tf remains more or less 
the same. Before this time, the change in Tf — T is due mainly to the variation of T. After t±, 
Tf begins to change very quickly. At t% the system has almost reached the structural equilibrium 
state. According to this behavior, the temperatures at those two times should correspond to the 
static and the dynamic T g . T(t\) = 103 °C should, therefore, be roughly the static T g . DSC 
results mentioned earlier suggest that this estimation is a few degrees beneath the actual value. 
Tfa) = 115 °C is approximately equal to the dynamic T g obtained by DSC. It should be taken 
into account that T g is not an input of the model. The only way that the model has to obtain the 
structural kinetics of the system is from the I data. For this reason parameters resulting from a 
fit to an experiment with T p = T po have been employed, since these parameters ought to be more 
representative of the behavior of the a relaxation as a whole. Nevertheless, it is remarkable that 
the picture of the glass transition is so realistic as the one that has been obtained. 

Now that we have seen that the results from the TNM model give a reasonable structural 
kinetics of the system, we are able to discuss which features characteristic of the TNM model or 
incompatible with it can be found in data coming from usual TSDC experiments. An important 
feature of the model is that it takes into account the structural state of the sample. For this 
reason, the results given by the model should depend on the thermal history. Nevertheless, the 
existence of a distribution of relaxation times can also yield this effect. In fact, a complete model 
should include some kind of relaxation time distribution modelization to be able to predict results 
coming from a wide range of experiments. Nevertheless, there are several ways to implement such 



modelization and a complete discussion lies out of the scope of this work. Instead, we will discuss 
several experiments and explain why their results are compatible with a TNM relaxation time 
even when the TNM model by itself may not be able to explain all the features observed in the 
experimental data, without making further assumptions. 

A usual way to obtain thermal history dependent results is changing the sub T g annealing time. 
The TNM model predicts that the structural state of the system at the end of a cooling ramp 
will depend almost exclusively on the cooling rate. As a consequence, the structural state of the 
sample during the TSDC ramp will vary in the same way in all the experiments. No change should 
be expected on the ln(r) plot due to this cause. This does not mean that such changes are not 
possible. In fact, polarization takes place when the sample is at different structural states. In 
the same way that the total polarization of the sample decreases with annealing time, the part of 
the mechanism that is polarized can also change. In fig. |H1 we present experimental results using 
annealing times between 3 min and 243 min. Prom these experimental data we can establish that 
in the case of PMMA the ln(r) plot does not change significantly during the first hour of annealing 
at 95 °. We can infer that we are polarizing the same part of the mechanism using this range of 
annealing times. The ln(r) plots change slightly for larger annealing times due to a small shift 
towards greater temperatures of the / plot. This does not mean that the TNM model is no longer 
valid but rather that we are not polarizing the same part of the mechanism. 

From the previous lines follows that there are only three parameters that we can change in order 
to observe a characteristic TNM behavior in the experimental data: storage temperature, cooling 
rate and heating rate of the TSDC ramp. 

The easiest way to obtain different values of Tf — T at the beginning of the TSDC ramp is 
employing several storage temperatures. Nevertheless, the structural state of the system will not 
vary very much neither at the end of the cooling ramp, during the storage time nor at the beginning 
of the TSDC ramp. This happens because temperatures suitables for storage are so far from T g 
that the structure of the system is almost frozen. The predicted intensity plots according to the 
TNM model can be seen at fig-HJ The hypothetical parameters that have been employed have been 
chosen in order to make the variation of the plots more noticeable. The foot of the a peak tends to 
grow for higher storage temperatures while the peak maximum exhibits a slight shift towards lower 
temperatures. This behavior is at a first glance different from the one expected from the Arrhenius 
or the WLF models, that predict no change in the data. Anyway, the TNM model can give a 
behavior similar to the one of the equilibrium models. The key parameter is the structural pre- 
exponential factor. For values of t s q less than 10~ 19 s, the difference between the plots becomes 



very difficult to appreciate and smaller than the experimental resolution. Experimental I plots 
obtained with two storage temperatures are presented in fig, 1101 The experimental data presents no 
differences between the two thermograms but this result is compatible with the parameters resulting 
from the fit to / data. Fits to / data may be more representative because the performance of the 
model is worse for structural states further from equilibrium and fits to / data tend to overestimate 
the influence of the data points placed near T^nax ? that are not as far from structural equilibrium 
as the other ones. 

Calculations assuming t s q ~ 10 -19 s show that few changes should be expected using different 
cooling rates. In fact, experimental results (not reported in this paper) show a displacement of the 
a peak towards higher temperatures as slower cooling rates are employed. This effect is explained 
easily if the existence of a distribution of relaxation times is assumed. It is not incompatible with 
a TNM behavior but it will certainly mask any characteristic TNM effect that could be observed. 

Finally, the heating rate of the TSDC ramp has a strong and well-known influence on the / 
data plots. The shift to higher temperatures of the a peak when the heating rate is increased does 
not depend on the model and can be explained using eq. El exclusively. The ln(r) plot will show 
almost no change, according to all the models. 

IV. CONCLUSIONS 

In this work we have shown how the TNM model can be used to modelize a more realistic shape 
of the TSDC intensity curve. Although the Arrhenius and the WLF models are the currently 
employed models in the literature, there are reasons to think that their description of the dielectric 
a relaxation is incomplete. The main one is that the structural state of the system should be taken 
into account in the modelization of a relaxation that is the dielectric manifestation of the glass 
transition. 

An important question is which is the relationship between the structural relaxation time given 
by the TNM model and the dielectric relaxation time obtained by the analysis of TSDC data. Our 
calculations indicate that a proportionality relationship between these two quantities is a useful 
hypothesis. The numerical results obtained in this way are able modelize the a relaxation over a 
wider temperature range than the Arrhenius or the WLF models. 

The application of the TNM model to dielectric data has its rough edges. The TNM model is, 
in fact, an extension of the Arrhenius model and therefore it is not surprising that its performance 
is better when the structural state of the system is not too far from equilibrium. As a consequence, 



the data points obtained at lower temperatures are not so well modelized as the ones placed nearer 
to the T g and the result of the fits depend somewhat on which is the magnitude for which the 
maximum likeness is searched. 

Nevertheless there are also some promising features. The TNM model provides a good de- 
scription of the relaxation over a wider range of temperatures than the equilibrium models. The 
modelization of the dielectric data provides also a realistic kinetics of the structure of the sys- 
tem, without the need of stating explicitly the glass transition temperature. The dependence of 
the available experimental data on the thermal history can be given by the TNM model if the 
pre-exponential factor of the structural relaxation time has values below « 1CT 19 . The parameter 
values obtained fitting the TNM model with an additional parameter to intensity data is compatible 
with the behavior observed in available experimental results. 

All in all, the numerical results support the convenience of taking into account the out of 
equilibrium state of the system, although a really fine tuned experiment would be needed to find 
qualitative evidence of its influence. 

This research is being continued along several lines. Other phenomenological models can be 
compared to the TNM model and may provide better results. The combination of phenomenological 
models with a distribution of relaxation times should give an improved description of the system. 
The use of other techniques such as dielectric electrical analysis can also give useful information. 
Finally, these new methods can be applied to review the classical studies on TSDC and physical 
aging. 
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TABLE I: Parameters of the experiments that have been performed. 



Experiment # T a (°C) t a (s) T p (°C) t p (s) V p (V) T s (°C) t s (s) 



1 


95 


1080 


95 


900 


800 


40 


840 


2 


107.5 


180 


107.5 


900 


800 


30 


840 


3 


95 


180 


95 


900 


800 


40 


840 


4 


95 


1080 


95 


900 


800 


40 


840 


5 


95 


1980 


95 


900 


800 


40 


840 




95 


3780 


95 


900 


800 


40 


840 


7 


95 


7380 


95 


900 


800 


40 


840 


8 


95 


14580 


95 


900 


800 


40 


840 


9 


95 


1080 


95 


900 


800 


40 


840 


10 


95 


1080 


95 


900 


800 


60 


840 



TABLE II: Results obtained from the fits to the Arrhcnius model 



Fit Experiment # Variable Range (°C) E a (eV) r (s) 

(a) 1 / 95.5-110.5 1.96 1.03 x 10" 

(b) 1 ln(r) 95.5-110.5 1.69 4.84 x 10" 

(c) 1 I 70.0-110.5 1.90 6.21 x 10" 

(d) 1 ln(r) 70.0-110.5 1.37 1.13 x 10" 



TABLE III: Results obtained from the fits to the WLF model 



Fit Experiment # Variable Range (°C) E w (cV) t (s) (K) 

(e) I I 95.5-110.5 2.60 x 10~ 2 2.70 x 10" 1 330 

(f) 1 ln(r) 95.5-110.5 2.60 x I0~ 2 3.45 x 10" 1 329 



TABLE IV: Results obtained from the fits to the TNM model 



Fit Experiment # Variable Range (°C) E a (eV) x t s0 (s) t^o (s) 



(g) 


1 


I 


70.0- 


-110.5 


2.08 


0.855 2.05 x 10 


-26 


T~sO 




(h) 


1 


ln(r) 


70.0- 


-110.5 


2.03 


0.607 9.53 x 10 


-26 


T~sO 




(i) 


1 


I 


70.0 


110.5 


1.52 


0.609 2.18 x 10 


-20 


1.36 x 10" 


19 


(j) 


1 


ln(r) 


70.0- 


-110.5 


1.13 


0.391 3.48 x 10 


-15 


2.09 x 10~ 


■14 


(k) 


2 


/ 


97.5- 


-117.5 


1.91 


0.462 2.27 x 10 


-24 


2.36 x 10" 


23 



TABLE V: Parameters of the experiments that have been simulated. 



E a (cV) X t s0 (s) T d0 (s) 

1.12 0.388 2.48 x 1CT 14 1.49 x 1(T 13 
Simulation # T a (°C) t a (s) T p (°C) t p (s) V p (V) T s (°C) t s (s) 



1 


95 


1080 


95 


900 


800 


30 


840 


2 


95 


1080 


95 


900 


800 


40 


840 


3 


95 


1080 


95 


900 


800 


50 


840 


4 


95 


1080 


95 


900 


800 


60 


840 



Figure 1 Outline of TSDC experiments. During the polarization time 800 V are applied to the 
sample. All the experiments begin and end at 125 °C or above. 

Figure 2 TSDC scan of the sample and its Bucci-Fieschi-Guidi plot. The following temperatures 
are marked in the figure: T x = 70 °C, T[ = 95.5 °C, T max = 104 °C and T 2 = 110.5 °C. 

Figure 3 Arrhenius model fits. The model values obtained are listed in table [H] Fits performed 
between T[ and T 2 : point-dash-dashed plots are fitted to I data (fit a). Point-dashed plots are 
fitted to ln(r) data (fit b). Fits performed between T\ and T 2 : continuous line plots are fitted to 
/ data (fit c). Dashed line plots are fitted to ln(r) data (fit d). 

Figure 4 WLF model fits. The model values obtained are listed in table Mil Continuous line 
plots are fitted to / data (fit e). Dashed line plots are fitted to ln(r) data (fit f). Fits performed 
between T[ and T 2 . 

Figure 5 TNM model fits. The model values obtained are listed in table IIVI Continuous line 
plots are fitted to / data (fit g). Dashed line plots are fitted to ln(r) data (fit h). Fits performed 
between T\ and T 2 . 

Figure 6 TNM model fits with additional parameter. The model values obtained are listed in 
table llVl Continuous line plots are fitted to / data (fit i). Dashed line plots are fitted to ln(r) (fit 
j). Fits performed between T\ and T 2 . 

Figure 7 Continuous line: calculated deviation from structural equilibrium (Tf — T) during ex- 
periments with different structural states at the beginning of the cooling ramp (t = 0). From this 
point, the thermal history of is the same as in experiment #2. Dashed line: Temperature of the 
sample (T). The following points are marked in the figure: T(ti) = 103 °C and T(t 2 ) = 115 °C. 
The parameters employed are the same as obtained in fit k. 

Figure 8 ln(r) results obtained with different annealing times. The experiments are listed in 
table. □ #3 (t a = 180 s, small circle), #4 (t a = 1080 s, large circle), #5 (t a = 1980 s, plus) and #6 
(t a = 3780 s, cross). The dashed line is the mean of the four data plots and is included as a guide 
for the eye. The ln(r) data were calculated from the experimental / data presented in the small 
frame, that also includes data from experiments #7 and #8. Lower peaks correspond to larger 
annealing times. 



Figure 9 Dependence on storage temperature for the TNM model with additional parameter. 
The simulations are listed in table E #1 (T s = 30 °C, continuous), #2 (T s = 40 °C, dashed), #3 
(T s = 50 °C, dot-dashed-dashed) and #4 (T s = 60 °C, dot-dashed). 

Figure 10 Experimental TSDC thermograms where two different storage temperatures have been 
employed. The experiments are listed in table □ #9 (T s = 40 °C, continuous line) and #10 
(T s = 60 °C, circle) 
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